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ABSTRACT 

Several  models  have  been  developed  in  the  past  to  describe  the  distribution  of  hydrometeor  canting  angles 
and  the  resulting  effects  on  polarimetric  radar  measurements.  These  models  are  reviewed  to  compare  their 
characteristics  and  to  assess  their  validity  under  specified  conditions.  Of  particular  interest  is  the  comparison 
of  the  traditional  two-component  model  and  a  two-dimensional  Gaussian  distribution.  Calculations  based  on 
the  Gaussian  distribution  are  presented  and  used  to  deduce  shape  distribution  parameters  from  radar  measure- 
ments.  ^  /  u  ■  J  c**  •  ^  '  7  -  T7  t,'  ■  v  .  /  "  ,  ■  . 


3^ 
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1.  Introduction 

The  distribution  of  canting  angles  is  one  of  the  mi¬ 
crophysical  attributes  of  hydrometeors  that  affect  the 
polarization-dependent  parameters  measurable  by  ra¬ 
dar  and  that  may  therefore  be  deduced  from  polari¬ 
metric  radar  measurements.  Recent  work  on  this  sub¬ 
ject  has  included  some  inconsistencies,  unnecessary 
complications,  and  erroneous  interpretations.  The 
fundamental  issues  which  need  to  be  addressed  in  re¬ 
lation  to  hydrometeor  canting  angles  are:  (i)  specifi¬ 
cation  of  a  form  of  the  canting  angle  distribution,  both 
in  one  dimension  and  in  two  dimensions,  and  (ii)  de¬ 
termination  of  reasonable  values  of  the  distribution 
parameters  for  various  types  of  hydrometeors.  In  sec¬ 
tions  2  and  3  we  discuss  key  developments  in  the  mod¬ 
eling  of  canting  angle  distributions  and  describe  radar 
measurements  which  support  theoretical  results  per¬ 
taining  to  the  standard  deviation  of  canting  angles  of 
raindrops.  Concepts  of  a  two-dimensional  distribution 
of  canting  angles  are  discussed  in  section  4,  and  a  com¬ 
parison  of  the  traditional  two-component  model  and 
the  two-dimensional  Gaussian  model  is  presented  in 
section  5.  Calculations  based  on  the  two-dimensional 
Gaussian  model  are  used  in  section  6  to  derive  esti¬ 
mates  of  mean  backscatter  amplitude  and  power  ratios 
from  polarimetric  radar  measurements. 

The  angular  parameters  of  hydrometeor  canting  are 
illustrated  in  Fig.  1 .  We  use  6  to  denote  the  absolute 
canting  angle,  which  is  the  angle  between  the  local  ver¬ 
tical  and  the  hydrometeor  symmetry  axis,  independent 
of  azimuth.  Following  a  long  standing  tradition,  we 


the  projection  of  6  onto  a  plane  perpendicular  to  the 
line  of  sight.  For  convenience,  most  of  the  following  L 
discussion  is  based  on  the  assumption  that  5^0.  ,  ^ 

2.  Theory  for  the  one-dimensional  distribution  ►  -  ' 

* 

Beginning  with  the  work  of  Saunders  (1971)  there 
has  been  a  tradition  of  treating  the  canting  angle  dis-  ' 
tribution  as  one-dimensional  in  a  plane  perpendicular  . 
to  the  radar  line  of  sight.  McCormick  and  Hendry 
(1975)  and  others  have  adopted  a  two-component 
model  of  the  apparent  canting  angle  distribution  in  (  v 
which  one  fraction  of  the  hydrometeors  has  a  single 
fixed  apparent  canting  angle  and  the  complementary 
fraction  has  a  uniform  random  distribution.  This  model 
of  the  distribution  function  was  used,  for  example,  by  C, 
Metcalf  ( 1 986)  in  the  computation  of  polarization  de¬ 
pendent  parameters  measurable  in  rain.  McCormick 
and  Hendry  related  the  effective  fraction  of  oriented 
scatterers  p„  to  a  distribution  function  by  the  equation 

J+w/2 

COs2(a  —  a)p(a  —  a)d(a  —  a),  ( 1 )  , " 

•»/ 2 

where  p(a  -  a)  is  the  distribution  function  relative  to 
the  mean  apparent  canting  angle  a.  They  did  not, 
however,  specify  the  form  of  the  distribution. 

From  the  two-component  model  the  reflectivity- 
weighted  average  of  the  complex  backscatter  amplitude 
ratio  defined  for  circular  polarization  is  given 

by 

re*,±2a>  =  Paoe*it2°\  (2) 


use  «  to  denote  the  apparent  canting  angle,  which  is  (McCormick  and  Hendry  ,975)  where  ,  is  the 

scatter  differential  phase  shift,  p„  is  given  by  Eq.  (1), 

-  and  the  upper  and  lower  signs  denote  the  transmission 

Corresponding  author  address:  Dr.  James  Metcalf,  AFGL/LYR,  of  right  and  left  circularly  polarized  signals,  respec- 
Hanscom  Air  Force  Base,  MA  01731,  tively.  The  quantity  veii  is  approximately  real  for  Ray- 
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PERPENDICULAR  TO  LOS 
IN  VERTICAL  PLANE 


IN  HORIZONTAL  PLANE 

Fig.  1.  Sector  of  the  hemisphere  on  which  pg  is  defined,  showing 
the  radar  elevation  angle  <t>,  the  canting  angle  6  of  the  symmetry  axis 
of  a  particle  from  the  vertical,  the  projection  a  of  the  canting  angle 
onto  the  plane  perpendicular  to  the  radar  line  of  sight  (LOS),  and 
the  angle  y  of  the  symmetry  axis  from  this  plane. 


leigh-Gans  scatterers  and  is  positive  and  negative  for 
prolate  and  oblate  spheroids,  respectively.  The  effect 
of  radar  elevation  angle  on  v  is  included  through  the 
assumption  that  for  an  individual  hydrometeor  the  ap¬ 
parent  amplitude  ratio  is  given  by  va  =  v0  cos2y,  where 
v0  is  the  amplitude  ratio  of  the  hydrometeor  when 
viewed  perpendicularly  to  its  symmetry  axis  and  y  is 
the  angle  of  the  symmetry  axis  from  a  plane  perpen¬ 
dicular  to  the  radar  line  of  sight  (see  Fig.  1).  If  the 
average  orientation  of  the  symmetry  axes  is  vertical, 
then  y  is  equal  to  the  radar  elevation  angle  <t>  and  v 
=  So  cos20.  Because  the  function  cos2a  passes  through 
a  full  cycle  in  the  interval  — x/2  «  a  «  jt/2,  the  ran¬ 
domly  oriented  fraction  (1  —  p„)  contributes  nothing 
to  the  average  in  Eq.  (2),  and  the  effect  of  the  apparent 
amplitude  ratio  is  evident  only  in  the  elevation  angle 
dependence  of  S.  The  reflectivity-weighted  average  of 
the  backscatter  power  ratio  v2  is  given  by 


v2^ 


p„  cosV  + 


(n)“ " 


(3) 


where  the  factor  ®^5  results  from  integrating  the  function 
cos4y  over  a  hemisphere  of  solid  angle.  Equation  (3) 
implies  an  ambiguity  in  the  two-component  model,  as 
the  distribution  over  a  hemisphere  is  inconsistent  with 
the  idea  of  a  one-dimensional  canting  angle  in  the  plane 
perpendicular  to  the  line  of  sight.  Specifically,  the 
quantity  S  in  Eq.  (2)  is  an  average  overj  distribution 
of  actual  shapes,  whereas  the  quantity  v2  in  Eq.  (3)  is 
an  average  over  distributions  of  both  actual  and  ap¬ 
parent  shapes. 


Beard  and  Jameson  (1983)  showed,  on  the  basis  of 
turbulence  theory,  that  a  Gaussian  distribution  of  ve¬ 
locity  fluctuations  would  yield  a  Gaussian  distribution 
of  the  tangent  of  the  apparent  canting  angle,  from 
which  one  can  derive  p(a).  [There  is  an  error  in  their 
equation  defining  p(«),  in  that  the  argument  of  the 
exponential  function  should  be  -tan2a/(2tri„„),  using 
the  present  notation.]  Their  calculations  showed  that 
the  standard  deviation  oa  of  the  apparent  canting  angle 
is  less  than  4°  for  a  thunderstorm  raindrop  size  distri¬ 
bution  and  turbulent  energy  dissipation  rates  up  to  0.2 
m2  s-3.  For  small  standard  deviations  of  apparent  cant¬ 
ing  angle,  the  differences  between  a  Gaussian  distri¬ 
bution  of  a  and  a  Gaussian  distribution  of  tana  are 
slight.  However,  for  large  values  of  the  standard  devia¬ 
tion,  the  Gaussian  distribution  of  tana  yields  unrealistic 
results,  as  it  does  not  yield  a  uniform  distribution  as 
<run«  approaches  infinity  and  it  yields  zero  probability 
density  at  a  =  ±t/2  for  all  values  of  a —  It  is  therefore 
useful  to  consider  the  properties  of  a  Gaussian  distri¬ 
bution  of  the  apparent  canting  angle.  Hendry  et  al. 
(1976)  and  Torlaschi  et  al.  (1984)  used  such  a  distri¬ 
bution  but  truncated  it  at  ±x/2  to  reflect  the  fact  that 
the  distribution  is  physically  limited  to  this  interval. 
The  resulting  relationship  of  oa  to  p„  is 

pa  =  exp(-2ffa2)  Re{erf[x/(2V20  +  /V2a„]}  (4) 

where  a„  is  expressed  in  radians.  Metcalf  and  Ussailis 
(1984)  suggested  that  the  distribution  in  the  interval 
±x/2  could  be  viewed  as  a  summation  of  segments  of 
an  infinite  Gaussian  distribution  folded  back  on  itself 
within  this  interval,  i.e., 

p(a)  =  (V&O'1  Z  exp[-(a  +  nx)2/( 2<ra2)\.  (5) 

CC 

Because  the  function  cos2a  is  cyclical  in  this  interval, 
the  limits  of  integration  in  Eq.  ( 1 )  can  be  extended  to 
infinity,  with  the  result  that 

p„  =  exp(— 2  <ra2)  (6) 

which  is  the  asymptotic  limit  of  Eq.  (4)  for  small  aa. 
The  approach  of  Metcalf  and  Ussailis  also  yields  the 
aesthetically  satisfying  result  that  dp{a)/da  =  0  at  a 
=  ±x/2. 

It  should  be  noted  that  neither  in  the  approach  of 
Hendry  et  al.  and  Torlaschi  et  al.,  nor  in  the  approach 
of  Metcalf  and  Ussailis  is  the  parameter  <ra  of  the 
Gaussian  distribution  function  exactly  equal  to  the 
standard  deviation  of  the  apparent  canting  angle  dis¬ 
tribution  defined  on  the  finite  interval  from  r/2  to 
+x/2.  The  latter,  denoted  by  aa,  approaches  the 
former  at  small  values  of  oa,  but  as  oa  approaches  in¬ 
finity,  the  resulting  uniform  distribution  on  the  finite 
interval  yields  a  standard  deviation  of  */V\2  radians, 
or  about  52°.  Hence  it  may  be  preferable  to  use  p„ 
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rather  than  a„  as  the  primary  descriptor  of  the  orien¬ 
tation  state  of  a  backscatter  medium. 

3.  Measurements 

Hendry  et  al.  (1987)  developed  a  method  of  deriving 
a  canting  parameter  from  measurements  with  varying 
linear  polarization.  They  showed  that  measurement  of 
the  maximum  and  minimum  linear  cross-polarized 
power  as  the  linear  polarization  vector  is  rotated 
through  at  least  180°  yields  the  parameter 

J+*-/2 

COs4(a  —  a)p(a  —  a)d(a  —  a),  (7) 

■w/2 

defined  on  the  basis  of  a  one-dimensional  distribution 
of  apparent  canting  angle,  which  differs  from  p„,  de¬ 
fined  by  Eq.  ( 1 ),  only  by  the  coefficient  of  the  argument 
of  the  cosine.  The  relation  of  p4  to  p„  depends  on  the 
distribution  function.  The  two-component  model 
yields  p4  =  p„,  and  a  one-dimensional  Gaussian  dis¬ 
tribution  function  yields  p4  =  p*.  Two  other  distri¬ 
bution  functions,  triangular  and  rectangular,  which 
seem  less  realistic  than  the  Gaussian,  yield  relations  of 
p4  and  pa  that  differ  only  slightly  from  that  of  the  Gaus¬ 
sian.  The  distinction  between  p4  as  a  measurable  pa¬ 
rameter  of  the  canting  angle  distribution  and  p„  as  a 
conceptual  parameter  highlights  the  importance  of 
specifying  a  physically  realistic  canting  angle  distri¬ 
bution  for  the  purpose  of  interpreting  radar  measure¬ 
ments.  Measurements  by  Hendry  et  al.  in  rain  yielded 
p4  =  0.914,  which  implies  <r„  =  6.1°  for  a  Gaussian 
distribution.  This  value  is  comparable  to  the  theoretical 
value  deduced  by  Beard  and  Jameson. 

4.  Theory  for  the  two-dimensional  distribution 

The  effect  of  canting  out  of  the  plane  perpendicular 
to  the  radar  line  of  sight  has  until  recently  been  included 
implicitly  with  the  effect  of  a  distribution  of  hydro- 
meteor  shapes,  as  might  result  from  raindrop  oscillation 
or  from  the  presence  of  multiple  types  of  hydrometeors. 
Recent  efforts  (Jameson  1987)  have  been  undertaken 
to  analyze  two-dimensional  canting  relative  to  the  local 
vertical,  with  the  goals  of  (1)  correctly  separating  the 
effects  of  shape  distribution  from  the  effects  of  canting 
angle  distribution  and  (2)  explicitly  describing  the  ap¬ 
parent  changes  of  canting  and  shape  distributions  in 
radar  observations  at  nonzero  elevation  angles.  The 
two-dimensional  distribution,  defined  on  a  hemi¬ 
spherical  surface,  must  be  represented  by  a  probability 
density  per  unit  solid  angle,  here  denoted  by  pa-  The 
distribution  is  assumed  to  have  a  maximum  value  at 
a  single  direction,  typically  near  the  vertical,  and  to  be 
axisymmetric  about  the  direction  of  the  maximum.  If 
the  maximum  of  pa  is  at  the  vertical,  then  p0  can  be 
expressed  as  a  function  of  8  only.  If  the  maximum  of 
Pa  is  displaced  from  the  vertical,  as  might  occur  as  a 
result  of  wind  shear  or  electric  field,  then  the  distri¬ 


bution  is  not  a  function  of  6  alone  but  can  be  expressed 
as  a  function  of  a  subsidiary  coordinate  relative  to  the 
location  of  the  maximum.  For  convenience  in  the 
computations  to  be  discussed  in  section  5,  we  assume 
the  maximum  of  pa  to  be  at  the  vertical.  The  angular 
coordinates  pertinent  to  this  distribution  are  illustrated 
in  Fig.  1 .  Here  the  symmetry  axis  of  a  hydrometeor 
lies  at  an  angle  6  from  the  vertical,  projects  an  angle  a 
on  the  plane  perpendicular  to  the  line  of  sight,  and  lies 
at  an  angle  y  from  that  plane.  If  a  =  0,  the  effects  of 
nonzero  value  of  y  are  equivalent  to  the  effects  of  a 
nonzero  radar  elevation  angle  </>.  A  nonzero  mean  ap¬ 
parent  canting  angle,  i.e.,  a  +  0,  affects  the  angular 
argument  of  certain  radar-measurable  parameters  but 
does  not  affect  the  magnitudes  of  the  quantities  dis¬ 
cussed  below. 

The  two-dimensional  distribution  of  6  yields  the  one¬ 
dimensional  distribution  of  a  by  the  relation 

J+*/2 

PaW«,  7))  cosyrfy  (8) 

-i/2 

where  the  angles  are  related  by 

cos0  =  cosa  cos0  cosy  +  sin0  siny.  (9) 

A  two-dimensional  Gaussian  distribution 

Pa(8)  =  (2Wr'  exp[— 02/(2<re2)]  (10) 

is  particularly  advantageous,  as  it  yields  a  Gaussian 
distribution  of  a  with  <r„  =  cre-  Jameson  (1987)  used  a 
two-dimensional  Gaussian  distribution  in  calculating 
radar-measurable  parameters,  but  the  utility  of  his  re¬ 
sults  is  diminished  by  his  arbitrary  selection  of  a,  -  18°, 
which  is  much  larger  than  the  values  deduced  by  Beard 
and  Jameson  (1983)  or  implied  by  the  measurements 
of  Hendry  et  al.  (1987).  Jameson  also  presented  results 
based  on  a  one-dimensional  Gaussian  distribution  of 
6.  Aside  from  the  conceptually  dubious  merit  of  ap¬ 
plying  a  one-dimensional  distribution  function  to  a 
two-dimensionally  distributed  variable,  the  practical 
value  of  such  a  distribution  function  is  uncertain  at 
best,  as  it  yields  a  singularity  of  the  probability  density 
per  unit  solid  angle  at  6  =  0,  does  not  yield  spherical 
uniformity  as  a,  -*■  oo,  and  does  not  yield  simple  re¬ 
lationships  between  pj.6)  and  p(a)  or  among  <re,  a„ 
and  pa. 

5.  Comparison  of  canting  angle  models 

a.  Formulation 

Calculations  have  been  performed  with  the  two-di¬ 
mensional  Gaussian  distribution  to  compare  this  model 
with  the  two-component  model  described  in  section  2 
and  to  determine  the  relationship  of  pdd)  to  p(a),  the 
effects  of  truncating  the  distribution,  the  mean  values 
of  backscatter  amplitude  and  power  ratios,  the  circular 
polarization  cross  correlation  measurable  by  radar,  and 
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the  correct  normalization  of  the  Gaussian  distribution 
defined  on  a  spherical  surface.  The  distribution  within 
the  domain  0  <  $  <  ir/2  was  approximated  successively 
by  Eq.  (10),  by 

pu(tf)  =  (2x<r»2)~ 1  (cxpl—02 /(2<r,2)) 

+  exp{-(tf-x  J*/(W)1}  00 

and  by 

pjfi)  =  (2 tc,2)'1  2  cxp{-(9  +  rt»)J/(W)].  02) 

*•-2 

Within  the  finite  domain,  Eq.  (10)  represents  a  Gaus¬ 
sian  distribution  truncated  at  0  =  x/2.  Equation  (II) 
represents  a  Gaussian  distribution  truncated  at  6  =  x 
and  folded  back  on  itself  at  6  =  x/2.  Equation  (12) 
represents  a  Gaussian  distribution  truncated  at  6  =  2x 
and  quadruply  folded  back  on  itself  in  the  interval  0 
«£  0  <  x/2.  The  successive  approximations  were  used 
to  assess  the  effects  of  truncation  and  to  assure  accurate 
results  as  a,  increased. 

We  assume  that  the  apparent  amplitude  ratio  is  given 
by  va (7)  =  v0  cos2 >,  as  in  section  2.  For  comparison  of 
the  two  distribution  models,  it  is  convenient  to  define 
the  amplitude  ratio  factor 

fA(a,)  =  vjv0  =  cos2y  cos2a  (13) 

and  the  power  ratio  factor 

fAat)  =  va2/v02  ~  COS47.  (14) 

For  comparison  with  the  measurements  of  Hendry  et 
al.,  we  use  their  formulation  to  define  the  quantity 

P*a  =  l(M)  cos4a]/(va2/v02) 

=  cos4y  cos4o/cos47  (15) 

which  is,  for  a  two-dimensional  canting  angle  distri¬ 
bution,  the  analog  of  p«,  defined  in  Eq.  (7).  The  overbar 
in  Eqs.  (13),  (14),  and  (15)  denotes  averaging  over  a 
hemisphere  weighted  by  the  distribution  function  pa- 
The  form  of  p«g  differs  from  that  of  p«  because  the 
distribution  function  pa(0)  jointly  affects  both  the  ap¬ 
parent  canting,  which  is  specified  by  a,  and  the  appar¬ 
ent  shape,  which  is  dependent  on  7.  Hendry  et  al. 
(1987),  in  developing  the  parameter  p«,  separated  the 
average  of  v2  from  the  average  of  cos4a.  Using  a  two- 
dimensional  distribution  function,  with  the  traditional 
assumption  that  the  shape  and  orientation  distributions 
are  independent,  we  factor  the  parameter  v2  into  two 
averages,  one  over  the  distribution  of  actual  shapes  and 
one  over  the  distribution  of  apparent  shapes.  In  this 
case  it  is  the  average  of  v2  over  the  distribution  of  actual 
shapes  which  is  canceled  out  of  the  ratio  of  maximum 
to  minimum  linear  cross  polarized  received  power,  and 
the  measurements  can  be  shown  to  yield  p*c,  defined 
in  Eq.  (15),  rather  than  p«,  defined  in  Eq.  (7). 

In  a  backscatter  medium  comprising  particles  of  a 


common  shape,  the  two-component  model  yields  the 
amplitude  ratio  factor 

fA  =  pa  cos2tf  (16) 

and  the  power  ratio  factor 

fp  =  Pa  cos*<t>  +  ^Vl  -  p,)  (17) 

corresponding  to  Eq.  (3). 
b.  Calculations 

The  results  corresponding  to  a  radar  elevation  angle 
of  zero  and  particles  of  a  common  shape  are  summa¬ 
rized  in  Table  1.  The  calculations  show  that  at  small 
to  moderate  values  of  ae,  the  two-dimensional  Gaussian 
distribution  of  6  yields  a  nearly  Gaussian  distribution 
of  a,  with  a„  as  a,.  As  ae  increases,  o„  approaches  its 
asymptotic  limit,  as  discussed  in  section  2.  Truncation 
effects,  not  documented  in  Table  1 ,  are  evident  only  if 
the  truncation  point  is  less  than  about  3.5a#.  In  other 
words,  Eq.  (10)  is  a  valid  approximation  for  ae  25° 
and  Eq.  ( 1 1 )  is  valid  for  <je «  50°.  Equation  ( 1 2)  should 
be  valid  for  a,  «£  100°,  at  which  value  p{a)  varies  by 
only  0.6%  between  0  and  x/2.  The  selection  of  one  of 
these  three  equations  to  approximate  the  Gaussian 
model  in  the  present  calculations  is  dictated  by  the 
magnitude  of  <r*.  All  yield  identical  results  if  a, «  25°; 
Eqs.  (11)  and  (12)  are  equally  usable  if  a,  50°. 

At  zero  elevation  angle  the  amplitude  and  power 
ratio  factors  calculated  from  the  Gaussian  model  are 
smaller  than  the  respective  factors  calculated  from  the 
two-component  model.  (Note  that  the  parameter  pa  is 
the  amplitude  ratio  factor  for  the  two-component 
model.)  This  results  from  integrating  over  the  hemi¬ 
sphere,  i.e.,  both  in  a  and  in  7,  and  reflects  the  decrease 
of  the  functions  cos27  and  cos47  as  7  deviates  from 
zero.  The  calculated  amplitude  and  power  ratio  factors 
for  several  elevation  angles  are  shown  in  Figs.  2  and  3. 
The  amplitude  ratio  factor  based  on  the  Gaussian 
model  is  well  fit  by  the  function 

L  -  exp(-3<T„2)  cos2(0)  ( 1 8) 

for  o, «  40° .  The  power  ratio  factor  at  <j>  =  0°  is  fairly 
well  fit  by  the  function 

fP  =  exp(-4<r„2)  +  |-~:j[l  -  exp(-4<r«2)].  (19) 

One  is  tempted  to  model  the  power  ratio  factor  by 

fP  =  exp(-4<7„2)  cosV  +  ^ j[  1  -  exp(-4<r#2)].  (20) 

However,  for  <t>  ?  20°  this  formula  is  a  less  accurate 
approximation  to  the  calculated  values  than  is  Eq.  ( 1 7). 
Therefore,  for  purposes  of  interpreting  radar  measure¬ 
ments  it  may  be  preferable  to  use  a  value  of  the  power 
ratio  factor  calculated  from  the  Gaussian  model  rather 
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Table  1.  Integrated  parameters  calculated  from  two-dimensional  Gaussian  distribution  of  absolute  canting  angle.  Calculations  pertain  to 
spheroidal  scatterers  with  mean  vertical  orientation  of  symmetry  axes,  viewed  at  horizontal  incidence  (zero  elevation  angle). 


00 

(deg) 

0m 

(deg) 

P« 

*J» 0 
(Gauss) 

r.W 

Pc 

Gauss 

2-component 

Gauss 

2-component 

5 

5.006 

0.985 

0.977 

0.985 

0.993 

0.985 

0.988 

10.05 

0.940 

0.912 

0.944 

0.972 

0.939 

0.954 

15 

15.18 

0.869 

0.816 

0.887 

0.939 

0.867 

0.897 

20 

20.46 

0.775 

0.700 

0.825 

0.895 

0.771 

0.819 

30 

31.44 

0.544 

0.457 

0.713 

0.787 

0.542 

0.613 

40 

40.65 

0.323 

0.261 

0.633 

0.684 

0.328 

0.390 

50 

46.34 

0.169 

0.135 

0.585 

0.612 

0.177 

0.216 

60 

49.36 

0.080 

0.064 

0.550 

0.571 

0.086 

0.105 

80 

51.53 

0.014 

0.011 

0.531 

0.540 

0.015 

0.010 

than  to  use  the  analytical  approximations  of  Eqs.  (19) 
and  (20). 

For  backscatter  from  particles  of  a  common  shape, 
the  two-component  model  yields  a  circular  polarization 
cross  correlation  given  by 

Pc  =  Pa  COS20  j\pa  COS *<j>  +  ^yjj(l  -  p„)j  .  (21) 

The  cross  correlation  from  the  Gaussian  model  is 
the  ratio  of  the  integrated  parameters  vj»0  and 
("a2/" o2)l/2-  As  shown  in  Table  1 ,  it  is  consistently  lower 
than  that  derived  from  the  two-component  model,  ex¬ 
cept  at  a,  -  80°. 

In  general,  the  integral  of  p^d),  defined  in  Eqs.  (10), 
(11),  and  (12),  over  a  hemisphere  should  yield  a  value 
less  than  unity.  This  is  due  to  the  fact  that  for  a  given 
value  of  0,  the  incremental  area  on  the  hemisphere, 


Ojtdag) 


Fig.  2.  Amplitude  ratio  factor  vjv 0  as  a  function  of  standard  de¬ 
viation  of  canting  angles  with  radar  elevation  angle  as  a  parameter. 
Solid  and  broken  lines  denote  calculations  from  the  Gaussian  and 
two-component  models,  respectively. 


given  by  sin 8d6d{,  is  less  than  the  incremental  area  on 
a  plane,  given  by  6d6d£,  where  0  is  a  colatitudinal  or 
radial  coordinate  and  f  is  a  longitudinal  or  azimuthal 
coordinate.  As  <re  increases,  the  distribution  is  spread 
to  regions  of  larger  0,  where  the  area  effect  is  more 
significant,  and  we  find  that  the  integral  yields  a  value 
approximately  equal  to  exp(-0.30o»2),  with  a,  in  ra¬ 
dians.  All  of  the  quantities  derived  from  the  Gaussian 
distribution  and  shown  in  Table  1  or  in  Figs.  2  and  3 
are  normalized  by  the  magnitude  of  the  integral  of  the 
probability  density  function. 

c.  Prolate  spheroids 

The  foregoing  theory  and  calculations  apply  to 
spheroidal  hydrometeors  of  any  ellipticity,  provided 


CTg(deg) 


FK5.  3.  Power  ratio  factor  as  a  function  of  standard  devia¬ 
tion  of  canting  angles  with  radar  elevation  angle  0  as  a  parameter. 
Solid  and  broken  lines  denote  calculations  from  the  Gaussian  and 
two-component  models,  respectively.  At  an  etevation  angle  of  31.3° 
in  the  two-component  model  the  power  ratio  factor  is  independent 
of  o9.  In  the  Gaussian  model  the  elevation  angle  at  which  this  factor 
equals  its  asymptotic  value  ’/is  varies  slightly  with  a,. 
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that  they  have  a  single  most  probable  orientation  angle. 
Prolate  spheroids  under  the  sole  influence  of  aerody¬ 
namic  forces  tend  to  be  oriented  with  their  symmetry 
axes  horizontal  but  uniformly  distributed  in  azimuth. 
The  appropriate  distribution  in  this  case  is  one  having 
a  minimum  at  6  =  0  and  a  maximum  at  6  =  x/2,  i.e., 
at  the  edge  of  the  hemisphere.  The  resulting  value  of 
fA  is  negative,  due  to  the  dominant  contribution  of  the 
cos2«  factor  in  Eq.  (13)  when  a  is  near  x/2.  Because 
the  complex  backscatter  amplitude  ratio  ve*i±la)  is 
positive  for  vertically  oriented  prolate  spheroids,  the 
average  of  this  ratio  over  a  distribution  of  nearly  hor¬ 
izontal  orientations  is  negative.  This  is  consistent  with 
the  fact  that  horizontally  oriented  prolate  spheroids 
appear  similar  to  vertically  oriented  oblate  spheroids, 
l  for  which  fA  is  positive  and  veii  is  negative.  The  distri¬ 

bution  of  prolate  spheroids  was  modeled  in  the  interval 
0  <  6  <  x/2  by  the  function 

PdS)  =  (2x<r,2)-‘  exp[-(0  -  x/2)2/(2<t42)].  (22) 

Calculations  corresponding  to  zero  radar  elevation  an¬ 
gle  are  illustrated  in  Table  2.  In  this  case  the  asymptotic 
value  of  \fA\  is  0.5  as  a,  -*■  0  and  that  of fP  is  0.375; 
these  are  the  averages  of  cos2y  and  cos^y,  respectively. 
For  these  calculations  the  quantity  oa  is  redefined  to 
represent  the  standard  deviation  of  the  apparent  ver¬ 
tically  oriented  oblate  spheroids.  The  parameter  pa  is 
similarly  redefined.  The  relation  of  pa  and  aa  is  similar 
to  that  in  the  case  of  oblate  spheroids.  Because  of  the 
broad  azimuthal  distribution  in  the  present  case,  ba  is 
much  larger  than  a,  when  a,  is  small.  As  a*  increases, 
<t„  approaches  the  same  limit  as  for  oblate  spheroids. 
The  amplitude  ratio  factor  fA  is  significantly  less  than 
P„  in  this  case  because  of  the  wide  distribution  of  sym¬ 
metry  axes  in  the  y  coordinate.  The  power  ratio  factor 
fp  increases  with  ae  to  approach  its  asymptotic  value 
of 0.533.  Application  of  the  two-component  model  to 
horizontally  oriented  prolate  spheroids  yields  the  am¬ 
plitude  ratio  factor 

U  =  -^jp»  COS 2(t>  (23) 


and  the  power  ratio  factor 

fp  =  (f  jpa  COS *4>  +  ^j(l  -  Pa)  (24) 

which  yield  the  cross  correlation  shown  in  Table  2. 
The  respective  factors  lh  and  J/»  result  from  the  fact 
that  in  this  case  the  “oriented  fraction”  is  distributed 
uniformly  in  azimuth. 

6.  Application  of  the  two-dimensional  model 

The  quantities  fA  and  //.enable  the  estimation  of  the 
reflectivity-weighted  mean  and  mean  square  of  the  ac¬ 
tual  amplitude  ratio  of  a  distribution  of  hydrometeors 
from  measurements  of  the  cross  covariance  amplitude 
ratio  (CCAR  ■  W/W2  in  the  notation  of  McCormick 
and  Hendry  1975)  and  the  circular  depolarization  ratio 
(CDR),  respectively,  if  a  parameter  of  the  canting  angle 
distribution  is  measurable  separately.  For  Rayieigh- 
Gans  scatterers  the  mean  actual  amplitude  ratio  is  given 
by 

I  Vo  I  =  |  CCAR  I//,  =  |  ORTT  |  (CDR)I/2//,  (25) 

where  ORTT,  in  the  notation  of  Hendry  et  al.  (1976, 
1987),  is  the  measured  cross  correlation  between  the 
two  received  signals.  The  mean  square  actual  amplitude 
ratio  is  given  by 

v?  =  CDR///.  (26) 

where  the  elevation  angle  dependence  of  the  measured 
parameters  is  incorporated  in  the  values  of  fA  and  //>. 
Equations  (25)  and  (26)  incorporate  the  traditional  as¬ 
sumption  that  the  shape  and  orientation  distributions 
are  independent.  As  used  here,  fA  and  fP  are  weighted 
only  by  the  probability  density  p<; ,  assumed  to  be  size- 
independent,  while  v0  and  p02  are  weighted  by  the 
number  density  and  backscatter  cross  section  as  func¬ 
tions  of  hydrometeor  size.  The  variance  of  the  actual 
amplitude  ratio  is  a,2  =  v02  - 
The  results  shown  in  Table  3  are  based  on  the  pub¬ 
lished  data  of  Hendry  et  al.  (1987).  Their  parameter 


Table  2.  Integrated  parameters  calculated  from  two-dimensional  Gaussian  distribution  of  absolute  canting  angle.  Calculations  pertain  to 
prolate  spheroidal  scatterers  with  mean  horizontal  orientation  of  symmetry  axes,  viewed  at  horizontal  incidence  (zero  elevation  angle). 


v.V».2  _ Pc 


(deg) 

V* 

(deg) 

Pa 

*’«/<' 0 
(Gauss) 

Gauss 

2-component 

Gauss 

2-component 

i 

7.95 

0.972 

-0.500 

0.375 

0.379 

-0.816 

-0.789 

3 

13.78 

0.917 

-0.496 

0.376 

0.388 

-0.809 

-0.736 

5 

17.77 

0.862 

-0.489 

0.377 

0.397 

-0.796 

-0.684 

10 

25.02 

0.728 

-0.457 

0.383 

0.418 

-0.738 

-0.563 

15 

30.39 

0.605 

-0.410 

0.394 

0.438 

-0.654 

-0.457 

20 

34.66 

0.495 

-0.355 

0.408 

0.455 

-0.556 

-0.367 

30 

40.89 

0.325 

-0.249 

0.441 

0.482 

-0.376 

-0.234 

40 

44.69 

0.216 

-0.172 

0.468 

0.499 

-0.251 

-0.153 

50 

46.95 

0.150 

-0.121 

0.486 

0.510 

-0.174 

-0.105 

60 

48.31 

0.110 

-0.090 

0.498 

0.516 

-0.127 
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Table  3.  Hydrometeor  shape  distribution  parameters  derived  from  radar  measurements  of  Hendry  et  al.  (1987). 


Case 

\ 

2 

Precipitation  type 

3 

Heavy  rain 

6  Jul  84 

Melting  layer 

10  Nov  84 

Snow 

4  Mar  85 

Elevation  angle  (deg) 

4.7 

16.3 

7.0 

ft  (Measured) 

0.914 

0.607 

0.280 

jCCARl  (Measured) 

0.179 

0.0139 

0.0222 

CDR  (Measured) 

0.0422 

0.00158 

0.00631 

oblate 

prolate 

a,  (Present  model) 

6.0° 

13.7° 

22.4° 

20.0° 

fA 

0.961 

0.777 

0.631 

-0.350 

fr 

0.966 

0.793 

0.783 

0.413 

’<S 

-0.187 

-0.0179 

-0.0353 

0.0634 

>o2 

0.0436 

0.00200 

0.00806 

0.0153 

a. 

0.0941 

0.0410 

0.0826 

0.106 

p4,  reinterpreted  here  as  p«<3,  yields  a  value  of  a,  through 
the  two-dimensional  Gaussian  model.  Successive  cal¬ 
culations  with  different  values  of  aa  were  used  to  de¬ 
termine  the  particular  value  of  a,  that  yields  the  mea¬ 
sured  value  of  pag-  This  value  of  ae,  in  turn,  yields 
values  of fA  and  fF  to  be  used  in  Eqs.  (25)  and  (26). 
Case  1 ,  an  observation  of  heavy  rain,  yields  a  mean 
amplitude  ratio  equal  to  that  of  raindrops  of  0.73  axial 
ratio,  which  corresponds  to  4.6  mm  equivalent  spher¬ 
ical  diameter  with  a  Pruppacher-Pitter  size-shape  re¬ 
lation.  The  standard  deviation  a„  equal  to  about  half 
the  magnitude  of  the  mean  v0  in  this  case,  implies  that 
essentially  all  the  hydrometeors  are  of  oblate  shape, 
which  is  to  be  expected  in  rain.  Each  of  cases  2  and  3 
yields  a,  greater  than  the  magnitude  of  v0eiS.  Case  2  is 
an  observation  of  a  melting  layer,  and  the  data  used 
for  the  present  analysis  are  from  near  its  base,  about 
200  m  below  the  altitude  of  the  maximum  circular 
depolarization  ratio.  At  this  altitude  the  hydrometeors 
are  likely  to  be  at  least  partially  melted,  but  the  larger 
derived  value  of  at  and  smaller  values  of  other  param¬ 
eters,  compared  to  those  of  case  1,  imply  that  the  hy¬ 
drometeors  cannot  be  approximated  by  a  single  dis¬ 
tribution  of  oblate  spheroids.  Case  3,  an  observation 
of  snow,  yields  potentially  ambiguous  results  because 
of  a  nonzero  mean  canting  angle  and  because  of  the 
more  likely  role  of  prolate  spheroidal  scatterers.  The 
radar  measurements  yielded  a  =  10°,  and  the  com¬ 
ponent  of  canting  relative  to  the  plane  perpendicular 
to  the  line  of  sight,  although  unknown,  may  be  of  com¬ 
parable  magnitude.  Hence,  at  the  radar  elevation  angle 
of  7°,  y  may  lie  between  -3°  and  17°,  and  the  mea¬ 
sured  value  of  pag  corresponds  to  <re  between  21.8°  (at 
y  =  17°)  and  22.5°  (at  y  =  0°)  for  oblate  spheroids. 
The  tabulated  quantities  for  Case  3  are  based  on  the 


assumption  that  y  =  7°.  Calculations  for  oblate  and 
prolate  spheroids  are  shown  for  comparison.  The 
transposition  of  signs  of  the  quantities  fA  and  S0e'i  be¬ 
tween  the  two  interpretations  indicates  that  in  either 
case  the  average  horizontal  dimension  of  the  scatterers 
is  greater  than  the  average  vertical  dimension.  The 
standard  deviation  a.  exceeds  the  magnitude  of  v0e’1 
by  a  smaller  factor  in  the  prolate  case  than  in  the  oblate, 
suggesting  that  the  prolate  description  is  somewhat 
better  than  the  oblate,  although  neither  is  fully  satis¬ 
factory. 

For  either  oblate  or  prolate  hydrometeors,  one  can 
hypothesize  a  distribution  having  a  relatively  large 
number  of  nearly  spherical  shape  and  a  relatively  small 
number  of  highly  nonspherical  shape,  which  would 
yield  a  small  magnitude  of  v0elt  and  a  relatively  large 
value  of  a,.  Such  a  distribution,  requiring  that  all  the 
hydrometeors  be  either  oblate  or  prolate,  seems  less 
realistic  than  a  mixture  of  oblate  and  prolate  shapes, 
the  analysis  of  which  is  beyond  the  scope  of  this  paper. 
Measurements  of  cloud  physical  parameters  by  air¬ 
borne  instruments  are  essential  to  the  development  of 
realistic  shape  distributions  but  were  unavailable  in 
these  cases.  In  evaluating  cases  2  and  3,  one  must  realize 
that  the  real  part  of  the  amplitude  ratio  veKi±2a)  is  a 
bipolar  function,  being  positive  if  the  larger  apparent 
dimension  of  a  spheroid  is  vertical  and  negative  if  it  is 
horizontal.  Hence,  a  combination  of  shape  and  ori: 
entation  distributions  may  yield  a  magnitude  of  v0e‘i 
much  less  than  the  standard  deviation  a,.  It  is  likely 
that  the  shape  and  orientation  distributions  are  not 
independent,  contrary  to  the  traditional  assumption. 
It  is  possible  that  the  approximation  va  =  v0  cos2y  may 
not  be  generally  valid  for  ice  phase  hydrometeors.  Fi¬ 
nally,  the  role  of  propagation  effects,  unknown  in  these 
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cases,  must  be  considered  in  the  proper  interpretation  likely  to  be  correlated.  Future  efforts  must  address  this 
of  polarimetric  radar  measurements.  more  complicated  analytical  problem. 


7.  Conclusion 

The  derivation  of  canting  angle  parameters  from  ra¬ 
dar  measurements  is  an  essential  element  of  charac¬ 
terization  of  hydrometeor  types.  As  polarimetric  me¬ 
teorological  radar  techniques  become  more  compre¬ 
hensive  and  more  quantitative  it  is  desirable  that  the 
underlying  theory  become  more  realistic  while  main¬ 
taining  simplicity  of  form.  Because  of  its  advantages 
and  simplicity,  the  Gaussian  distribution  of  canting 
angles  (both  one-dimensional  and  two-dimensional) 
seems  to  be  a  desirable  standard  until  such  time  as 
there  is  a  strong  rationale  for  using  an  alternative  dis¬ 
tribution.  Theory  and  measurements  indicate  that  the 
standard  deviation  of  canting  angles  of  raindrops  is 
typically  about  4°-6°.  Differences  in  amplitude  and 
power  ratio  factors  derived  from  the  Gaussian  and  two- 
component  models  are  small  enough  in  some  cases 
that  they  do  not  substantially  affect  calculations  such 
as  those  of  Metcalf  ( 1 986).  On  the  other  hand,  the  dif¬ 
ferences  are  large  enough  to  justify  the  use  of  the  Gaus¬ 
sian  model  in  preference  to  the  two-component  model 
in  the  interpretation  of  radar  measurements. 

In  this  and  all  previous  work  on  this  subject,  authors 
have  assumed  the  independence  of  canting  angle  dis¬ 
tribution  and  particle  shape  distribution.  This  as¬ 
sumption  permits  separation  of  the  respective  physical 
parameters,  e.g.,  a,  p„  and  P0.  In  practice,  since  the 
distributions  of  canting  angle  and  of  shape  are  likely 
to  be  size-dependent,  the  respective  parameters  are 
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